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Abstract 

Application of fractional calculus to the description of anomalous diffusion and re- 
laxation processes in complex media provided one of the most impressive impulses 
to the development of statistical physics during the last decade. In particular the 
so-called fractional diffusion equation enabled one to capture the main features of 
anomalous diffusion. However the price for this achievement is rather high - the 
fractional diffusion coefficient becomes an involved function of a characteristic of 
the media (e.g., that of the radius of pores in the case of the porous one). Revealing 
this dependence from the first principles is one of the main problems in this field 
of science. Another one still remains that of extracting this dependence from the 
experiment. The latter problem is tackled in the present paper. Our aim is to pro- 
vide detailed and pedagogical deriving the relationship of the fractional diffusion 
coefficient with experimentally observable value from nuclear magnetic resonance 
(NMR) spin-lattice relaxation data. The result obtained promotes the NMR relax- 
ation method to become a powerful tool in solving the problem of experimental 
measuring the fractional diffusion coefficient. Also the merits and limitations of 
NMR relaxation method and pulsed-field gradient (PFG) NMR for the research of 
anomalous diffusion are compared and discussed. 
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1 Introduction 

Revealing the origin and adequate theoretical description of diffusion and re- 
laxation processes have been one of the main concerns for statistical physics 
from the very beginning of this field of science. The obstacles that prevent 
achieving these aims are augmented substantially with the increase of the 
complexity of the system under consideration. That is why the trend to re- 
search these processes in the so called complex systems (for review see PQ, [2], 
[3] and refs. therein) posed new challenges for statistical physics and gave new 
impetus to its development. The history of this trend goes back as far as to 
the middle of the 19-th century when Kolraush introduced his famous streched 
exponential function for the description of charge relaxation in a Leyden gas. 
However only the end of the 20-th century is marked by a real breakthrough 
and by an explosion of activity in this field of science. 

Fractional dynamics is a modern and fruitful approach to the description of 
anomalous diffusion processes in randomly disordered media (see [2], [3] and 
refs. therein). The idea of this approach goes back to many pioneers whose 
achievements are honored in the reviews mentioned above. The reason for 
its introducing is as follows. As is well known the mean squared displace- 
ment of a free particle in a homogeneous media grows linearly with time 
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< (r — < r >) 2 >= 2dC\t where C% is the diffusion coefficient with the di- 
mension cm 2 / s and d is the (embedding) spatial dimension. This conventional 
Einstein relationship of the classical theory results from the ordinary diffusion 
equation for the probability density function to find the particle at position 
x at time t and is a direct consequence of the Fick's second law. However 
this simple law does not take place in the inhomogeneous disordered media. 
The generalization of the Fick's law for complex systems requires conceptually 
new physical ideas and even envokes to new and rather unusual for common 
audience mathematics - the so called fractional calculus [3], [5], [B], PQ. 

One of the basic tools of the fractional dynamics is a fractional diffusion equa- 
tion (FDE) in which a fractional diffusion coefficient (FDC) C a with an un- 
usual dimension cm 2 / s a is presented. The case of subdiffusion (0 < a < 1) 
is ubiquitous in nature yielding the mean squared displacement < (r — < 
r >) 2 >= 2dC a t a and it originates in any fractal media due to the presence 
of dead ends on current ways. The FDC acquires a functional dependence on 
the phenomenological parameter a referring to the extent of inhomogeneity 
of the system with a = 1 corresponding to the case of ordinary diffusion. To 
establish this dependence for each class of systems is one of the main problems 
for the theory. 

On the other hand the problem has been posed to determine the FDC exper- 
imentally. It was pioneered by [7] where the interdiffusion of heavy and light 
water in a porous media was observed by means of NMR. The authors made 
use of a one- dimensional form of the FDE and found C a using the special case 
of its a = 2/3 solution. In [8], [9] the procedure is called "neither very accurate 
nor of general use". Though such critique may be justified for the particular 
procedure of the paper [7j one should not think that the abilities of NMR for 
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measuring FDC are limited somehow. From their own side the authors of [8], 
[9J proposed a method to measure the FDC with the help of a membrane sys- 
tem where the substance of interest is transported in a solvent from one vessel 
to another across a thin membrane. In our opinion namely their method is not 
of general use because it requires the incorporation of a membrane into the 
system of interest and thus can hardly be applied to, e.g., porous materials. In 
contrast NMR belongs to non-destructive and non-invasive methods and thus 
can indeed be of general use. NMR diffusometry [10] is a powerful method for 
investigation of subdiffusion processes and it is widely applied to exploring 
transport in porous [TO] , [TT] . [12], percolative [13], [14] . [7J and polymeric [T5] 
systems. The theory of NMR diffusometry in disordered media [11], [7J, [16] 
is developed within the framework of fractional dynamics [2], [3]. The aim 
of this chapter is to show how one can retrieve the functional dependence of 
the FDC on the parameter a from the NMR spin-lattice relaxation data. We 
obtain a simple analytical formula relating the FDC with the contribution to 
the spin-lattice relaxation time by anomalous translational diffusion that can 
in principle be extracted from the experiment. 



2 Setting the stage 

As is well known in the theory of spin-lattice relaxation by dipole-dipole in- 
teraction (see e.g. [1 TJ ) the contribution to the spin-lattice relaxation rate 
constant due to translational diffusion with the spectral density at a Larmor 
frequency u>l of the correlation function for spherical harmonics has the form 

(m) trans = 3 7W+i) {J( i) K) + j(2 ){2ul)} (1) 
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where 7 is the gyromagnetic ratio of the nucleus, I is their spin and h is the 
Planck constant. The spectral densities are proportional to the function J(ui) 
of the spectral density 



JW( W ) = 1 |JH (2) 



and 



j< a V) = |pM (3) 



Thus the function J(o>) becomes the key object of the analysis. The consid- 
eration of the spin-lattice relaxation in homogeneous media due to diffusion 
motion of the particles presented in [17] is based on the ordinary diffusion 
equation for the probability density. To describe the case of inhomogeneous 
and disordered (e.g., porous) media the generalization of the ordinary diffusion 
equation within the fractional calculus was suggested by Schneider and Wyss 
[T8] which includes the FDC C a with the dimension [cm 2 /s°]. We introduce 
the characteristic time 

where d is the least distance to which the molecules can approach to each 
other (if the molecules are considered as spheres of the radius a then d = 2a 
[T7]). For the homogeneous case the normal diffusion coefficient C\ is given by 
the Stokes formula 



1 Quarj ^ 
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where r\ the viscosity of the media and consequently the characteristic time is 



127ra 3 T? 



The generalization of the ordinary diffusion equation within the fractional 
calculus suggested by Schneider and Wyss [18] is 



^^ = C a (D^P)(r,t) (7) 



where P(r, t) is the probability density function to find the particle at position 
r at time t, V 2 is the three-dimensional Laplace operator, C a denotes the 
fractional diffusion constant with the dimension [cm /s a ] and Dl~ a is the 
Riemann-Liouville fractional derivative of order 1 — a and with lower limit 0+ 
which is defined via the following relationship [I] 

(^)w = f^s/ (I -» rl/to) * (8) 

v ' 



where is a gamma function. The solution of the equation (j7j) for the case 
of sub-diffusion < a < 1 with the initial condition P(r, 0) = <5(r) where S(x) 
is a Dirac function is obtained in [18J and expressed via the Fox's function 



P(r,t) 



)3/2 



H 



AC n t a 



'I, a) 



(3/2,1), (1,1) 



(9) 
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The latter is defined as H™ n 



\ 



(an, At), {a p ,A p ) 



j=lf L dsz s 7]{s) 



J 



(bt,Bt),...,(b p ,B p ) 

where n[s) = tip -n, — i a \ttS — ^r, i — J- he nomenclature in the r ox s 
function associated with the vertical bar is explained via its explicit definition 
by the contour integral. The requirements to the contour path L are formu- 
lated in EDI. 



3 Direct problem 



In this preliminary Sec. we solve the direct problem - how to calculate the spec- 
tral density J(to) knowing the FDC C a ? - and present the detailed derivation 
of the result obtained in [16]. We follow the algorithm of [17J (in what follows 
all corresponding results from [T7] are obtained as a particular case a=l of the 
present approach). Under r we denote the vector ri — T2 connecting two iden- 
tical molecules diffusing relative to each other rather than the radius- vector of 
the molecule diffusing relative to a fixed point. This leads only to the change 
of 4C a t a by 8C a t a in (J6]). Our aim is to calculate the correlation function 

G(1) = 4/5ip)5W|lip (t _ ri , 1)A/r (10) 

where N is the number of spins in 1 cm 3 , T™(0, (p) is a spherical harmonic, 
P(r,t) is given by (Q and * denotes complex conjugate. To be more precise 
we need the spectral density of the correlation function G(t) to calculate the 
spin-lattice relaxation rate constant with the help of (J2j). At integration in 
( 11 01) one should take into account that r and r can not be less than some 
limit value d - the least distance to which the molecules can approach to each 
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other. If the molecules are considered as spheres of the radius a then d = 2a 



First we make the Fourier transforming of the function P(r — r ,t) with the 
mentioned above change of 4C a t a by 8C a t a 



t 



r - r 



20 

3 -"12 



(r - r ) s 

8C a t a 



\ 



(l,a) 

(3/2,i),(i,r 



(2tt)< 



c? 3 v/(v) exp 



zv(r - r , 



Denoting 

R = r - r ; 



R = | r — r 



v = v 



(12) 



we have the inverse transform (making use of spherical coordinates) 



47T 

V 



/(«) = — / dR —^'"' h? 



sm(vR) 



R 2 



8C n t c 



\ 



\ 



'l,a) 



(3/2,i),(i,i; 



/ 



(13) 



Substituting f[T3"j) into ffTTj) . then substituting the result into flTUJ) and making 
use of the known identities for spherical functions 



da rf (n)Tf,(fi) = <w<W 

1/2 



exp(— zvrj 



47r (^:) E^Tf(n)Tf(fi')Jm/ 2 K) 



(14) 
(15) 



where J u (x) is a Bessel function of order v we obtain 

2 



G(t) 



3/2 



7T 



dv f(v) 



o 



Jb/2{vr) 



r 



3/2 



(16) 
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where f(v) is given by (fT3|) . Taking the known value of the inner integral 



dr^l = v^(vdr^J 3/2 (vd) 



(17) 



and denoting u = v d we obtain 



°fdu 




00 




J3 (u) 


7 


J 


La. 













JP Sm v^T ' X 



i? 2 



(l,a) 

(3/2,1), (1,1) 



;is) 



The integration can be fulfilled and yields the expressions for both the corre- 
lation function and its spectral density via Fox's functions. However a Fox's 
function is to regret not tabulated at present either in Mathematica or Maple 
or Matlab. Thus it is rather difficult to use such formulas for plotting the be- 
havior of the correlation function and the spectral density. That is why it is 
useful to obtain another representation of the spectral density which enables 
one to plot the frequency dependence of this function. For this purpose we 
make the Fourier transforming of the Fox function in fflBl . The latter can be 
achieved with the help of the following trick going back to original investiga- 
tions of Fox. First we make use of (35) from the paper [18] to write 



rr20 
H 12 



R 2 



8C a t c 



\ 



'l,cO 



(3/2,1), (1,1) 



1720 



Rt- a ' 2 
(8C^ 



(W2) 

(3/2, 1/2), (1,1/2) 



(19) 
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Then with the help of transform (57) <-> (56) from [19] we obtain 



:\i ( - if 



R 



-a/2 



(8C Q r«)V2 I Tq 



\ 



(l,a/2) 

(3/2, 1/2), (1,1/2) 



/ J 



2-2s/c 



2/? 



a l(8C Q r°)V2 



r( '3_£A r(l-s/a) 



2 a/ rri-s 



(20) 



where M{...; s} denotes the Mellin transform and r a is given by (j4j). We in- 
troduce the designations 



2R 



x 



(21) 



and make use of the identity [23] 



F c {/(x); ^ = M" 1 ] r(s) cos y M {/(*); 1 - s} 



(22) 



where i 7 ^ denotes the cosine Fourier transform and M denotes the in- 

verse Mellin transform. Thus for our function of interest 



/(*) = #12 



R 2 



8C a r' 



-X 



\ 



l,a) 



(3/2,1), (1,1) 



(23) 



we have 



M{f(x);l-s} 



2(l-s) 

1 fr\ a /3 1 sN l r fi 1 s 
T(s)a\2 I \2 a aj \ a a 



(24) 



and 



F {f{x); z} = 2F C {/(*); z) = ^ (0^ G{z) 



(25) 
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where 



W.r* „= (i p r (|_I + i) r(l _I + i)j (26) 



The calculation of the inverse Mellin transform in this formula requires routine 
manipulations with the help of its known properties [24J 
a~ s g(s) <r+ f(ax); g(s/h) ^ hf(x h ); T(s) <->• e~ x 

CO 

9l (s + A)g 2 (s + A + B + l)^x A Jd£ J\( r. Oh(0- 

o 

sin(avrs/2)r(s) <-> e -*™<™/2) s i n (x sin(avr/2)) ; 

cos(a7rs/2)r(s) <->• g~ x ' cos ( 7rQ / 2 ) qos (a; sin(ct7r/2)). As a result of tedious but 
straightforward calculations we obtain 

G(z) = ~a((£j 2/a z) x 



J f 7 exp(-0 exp ^ — — cos — j sin j — sin- — j (27) 



The substitution of the results into spectral density of the correlation function 
(HUD yields 



1 ; nVWC „(ut, n ) l ~ a J J 



J 3 (u) 

2 







2 

sin I — — ) x 



OO , 

J ^r 1/2 exp(-Oexpl- 

n V 



R 2 (ur n ) a 7ia\ 

cos — — X 



&C a rg 2 
. [i? 2 (u;r a ) a . na 7ia\ 

sin < sin > (2e 

1 f8C Q r£ 2 2 J 1 



We denote 



R 2 (uT a ) a na na R 2 {uiT a ) a ira 1 

a = ^T sm T ; 6 = T ; c= ^T cos T ; x = e (29) 



n 



Then the inner integral over £ is 

oo ^ 

I = cos b J dx x~ 3 ^ 2 exp ( — \ exp(— 
o 

oo _^ 

sin b J dx x _3//2 exp ( — \ exp(— cx) cos( 



cx) sin [ax — 



(30) 



Both integrals here can be calculated with the help of N2.5.37.2 from [25]. As 
a result of straightforward manipulations we obtain 



7i \ 1 / 2 ( 7ra\ ( ira na 

I = 2o — exp — h cos — sin [h sin 

' \2hJ V 4 J V 4 2 



(31) 



where 



1/4 



We denote 



B 



2R(uT a 



,a/2 



(32) 



(33) 



The integral over i? in flze 



/ ii \ / 71 Ot 

dR sin ( — — J exp I— BR cos — ^- ) sin ( BR sin 



7ra ira 



(34) 



after simple manipulations can be cast into the combination of the integrals 
of the type N2.5.30.8 from [25]. As a result we obtain after straightforward 
calculations 

du 3 sin 



K 



(35) 



■ft 



4C^ 



Taking into account the definition we finally obtain the formula for spectral 
density of the correlation function in the general case of inhomogeneus media 
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(0 < a < 1) PI] 




DC 



J(w) 



C a (i (wr tt ) 



7ra 
sin — 
2 



/ 



du J s (u ) x 







(36) 



4 + (^r a )' a + 2u 2 (u;r a ) a cosf 



This formula generalizes that VIII. 113 from [T7] referring to a particular case 
a = 1 to the case of arbitrary < a < 1. The integration in it can be easily 
done with the help of Mathematica or Maple or Matlab. Thus this formula is 
convenient for plotting the spectral density (see [16] ) - The substitution of this 
formula into (CQ) , (J2j) , ([3]) solves the direct problem, i.e., how one can calculate 
the contribution into spin-lattice relaxation time by anomalous translational 
diffusion (1/Ti) trans knowing the FDC C a l 

4 Inverse problem 

In this central Sec. we tackle some more difficult inverse problem, i.e., how one 
can calculate the FDC C a knowing the contribution into spin-lattice relaxation 
time by anomalous translational diffusion (l/Ti) tTOn , s ? Our solution of this 
problem is based on the fact that the integral in the (|28|) can be cast into a 
form of convolution for Mellin transform J °° ^(f(u)K «-> (p(s)K(s). We 
factorize the characteristic time into the dimensional part T\ ([ti] = s) and a 
dimensionless function f(a) (/(l) = 1 ) accounting for the dependence on the 
parameter a 



T a = Tif(a) 



(37) 



13 



The function f(a) plays a key role in our analysis because knowing it we have 
r a and from (j4]) we obtain the desired value of the FDC C a . We denote the 
dimensionless variable 

x = ujlTi (38) 



and introduce the function 

H^a) = -—^^——(^\ (,-.,.) (39) 



%nN^h 2 I(I + l)x VTj 



1/ t 



rans 



We consider this function as being extracted in principle from the experiment, 
i.e., as a primarily given one. The condition of normalizing for the function 
H(x, 1) by the requirement f(l) = 1 will be given below. Making use of the 
substitution u = v a ' 2 we obtain from (TTl) and (1361) 



2xH( x, a) = a(f(a)T sin — / — [j 3 (V /2 Y 2 

2 J f L 2 V / 


1 1 + (x/v) 2a (f(a)) 2a + 2 (x/v) a (f(a)) a cos f 
2 Q+1 (x/w) Q 

1 + (x/v) 2a (2f(a)) 2a + 2 (x/v) a (2f(a)) a cos f 



x 



(40) 



Now we apply Mellin transform in the variable x (which we denote M{f(x), s} = 
g(s) or <->) to both sides of this equation. Making use of the property x 13 f(ax h ) <- 
ft-i fl -(»+fl/^[(4+j3)/ft] where a > 0; /t > and of N6. 8. 33 from [21] we obtain 

[Ja (x^ 2 )] 2 - 2^ r(l-2 8 /a)r(3/2 + 8 / a ) 

J J a T 2 (1 -s/a)T (5/2 -s/a) 1 j 

where — 3a/2 < Re s < a/2. Making use of N6.2.12 from [21] we obtain 

x a 



1 + x 2a (f(a)Y a + 2x a (f(a)) a cos ■ 
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(f(a))- (s+a) vr sin(7rs/2) 



a sin(7ra/2) sin(7rs/a) 



(42) 



where —a < Re s < a. Thus both transforms (j4ip and (1421) have the common 
region 

-a<Res<a/2 (43) 



For ( 1421) to be valid we must be sure that 

(f(a)) a > (44) 

To verify the latter we obtain from (T4J) and ([3 
\d? ) \2C U/ 

We see that at C a be a decreasing function with the decrese of a from the value 
a = 1 (that is an inherent feature of the subdiffusion) the function (f(a)) a 
is an increasing one from the value (/(a = l)) a_1 = 1, i.e., indeed (144")) takes 
place. 

Applying (|4T!) and (142!) to the convolutions in (HUI) we obtain 

owr rr/ x i 2 2s /°r(l-2s/a)r(3/2 + s/a) 
2M{xi/(a;,a),s} - ' 



a r 2 (l-s/a)r(5/2-s/a) 

vr (/(a))~ s sin(7rs/2)(l + 2 1 - s ) 
sin(7rs/o;) 



(46) 



Making use of the properties of a T-function we can obtain after straightfor- 
ward calculations 

T 2 (1 - s/a) T (5/2 - s/a) _ tt2 4s / q - 1 (3/2 - s/a) (1/2 - s/a) 



T (1 - 2s/ a) T (3/2 + s/a) (i/ 2 + s /«) sin ( f ) V {2s/ 



'a) 
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Substituting it into (145!) we obtain 



sin (^) T (2a/ a) (2^f(a)] 

(3/2 - s/a) (1/2 - s/a) 



aM {xH(x, a); s} 



(1/2 + s/a) (1 + 2^) 



(4£ 



With the help of N7.3.9 from pi] we obtain 



sn. [ y ) T (2,/«) (2 2 /«/(a) 



a 

2 6XP 



-2(x /( a)r/ 2 cos(^) 



sin [ 2 (x/(«)) Q/2 sin ( ^ 



(49) 



where Re s > —a/2. Thus the common region for all transforms is 



a/2 < Re s < a/2 



(50) 



Taking the inverse Mellin transform for both sides of (|48p we obtain 



-exp 



-2(xf(a)) a/2 cos^ 



sm 



2(x/(a)r /2 sin^ 



duH(u, a)G ( — 



(51) 



where the fuction G(x) is defined as the inverse Mellin transform 
(3/2 - s/a) (1/2 - s/a) 



G(x) 



(l/2 + s/a)(l + 2 1 - 



(52) 



The poles of its righthand side are s = —a/2 and the solutions of the equation 
1 + 2 1 ~ s = that are given by s m = 1 - ^(2m + 1) where m — 0; ±1; ±2; .... 

Further we consider the case 



x < 1 



(53) 
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that is sufficient for all practical experimental situations. Then at calculation 
of the integral 



K J 2ni J (l/2 + s/a)(l + 2 1 " s ) 



(54) 



we can choose < o < 1/2 and close the contour in the lefthand halfplane 
(due to (1531 ). All the poles s m appear to be beyond its interior and only the 
pole s = —a/2 is inside it. Thus 

, (3/2 -s/a) (1/2 -s/a)l 2a;Q/2 



G(x) = res I x~ 



(l/2 + s/«)(l + 2 1 - s ) 



s=-a/2 



1 + 2 1+Q / 2 



(55) 



As a result we obtain 



exp 



-2 (xf(a)) a/z cos 



a/2 ™ 



sm 



2(xf(a)) a/2 sin 



7T« 



1 + 2l+°/2 



duH(u,a)u- a/2 



(56) 



This equation implicitly defines the required function /(a) via the given func- 
tion H(x,a). It is useful to cast it into another form making use of the fact 
that the left side of this equation has the form of the producing function for 
a Chebyshev polynom 



n=l 

4x a/2 
I + 2l+«/2 



duH(u, a)u 



-a/2 



(57) 



Now we take into account that in all practically important situations we can 
adopt the requirement 



x « 1 



(5f 
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(e.g., in [16] we had x = = 0.015). In this case we can restrict ourselves by 
the first term in the sum only. Taking into account that U\ (cos ™J = sin ^f- 
[26] we finally obtain 

2/a 



/(«) 



sin (?f) (1 + 2!W2) J 



duH(u, a)u 



-a/2 



(59) 



From the requirement /(l) = 1 we obtain the condition of normalizing for the 
function H(xA) 



J duH(u, IK 1 / 2 « 1 + 2- 3 / 2 
o 



(60) 



The desired dependence of the FDC on the parameter a hence has the form 

-2 



1-a 



sin (^f) (1 + 2 1 " 



-a/2 A 



duH(u, a)u 



-a/2 



(61) 



This formula is the central result of this chapter. 



5 NMR relaxation or PFG NMR? 



As is well known there are two strategies to employ NMR for measuring trans- 
lational diffusion coefficient: the NMR relaxation method [17] and the PFG 
NMR (for review see [27], [2B] and refs. therein). Both of them have their mer- 
its and limitations nicely discussed in [27]. As a matter of fact in the realm of 
ordinary diffusion the PFG NMR is a fine Prince while the NMR relaxation 
is a Beggar. However their roles are not so obvious for the case of anomalous 
diffusion. 

The theory of PFG NMR is developed on the base of a combination of Bloch 
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equations with an ordinary diffusion term pioneered by Torrey [29], [30]. The 
scheme works well for the case of unrestricted isotropic diffusion but becomes 
analytically intractable for the case of anisotropic diffusion or that within a 
confined geometry While simple geometries are still amenable to analytical 
approximations more complicated ones require numerical solutions [27]. It is 
rather problematic to apply this scheme to anomalous diffusion in complex 
media (in the sense of [I], [2], [3]) where even the geometry as a rule can not 
be defined. 



It is tempting to develop the theory of PFG NMR suitable for the description 
of the anomalous diffusion along the line of the papers [TT], [7J, [TB] and this 
chapter. However such a project also encounters severe difficulties. To ehxibit 
their origin we briefly scetch the mathematical scheme referring for explana- 
tion of designations to ,e.g., [30]. In our case instead of (G.7) from [30] we 
have a fractional Torrey equation 



dM+(r,t) 
dt 



-I'JZ 



'dH\ . M+(r,t) 
— I M + (r,t) 



dz I 



C a (Dl- a V 2 M + )(r,t) 



(62) 



Analogously (G.9) from [30] we seek the solution in the form 



M + (r. t) • ( »xp ( ~ ) exp 



• fdH\ ' 



A(t) 



(63) 



Then for the function A(t) we obtain a fractional differential equation 



oxp ! -— ] oxp 



• f9H\ ' 



dA(t) 
dt 



-^(7^) 2 ^{t 2 exp(-^)exp 



—i^z 



'dlf 

dz 



A(t) 



(64) 
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Recalling the definition of the fractional derivative © we can rewrite AMI) as 
an integro-differential equation 



OXj) [ -— J OX] ) 



- %1Z Ur r 



dA(*) n i dH \ 

— = - c « 7 57 x 



1 d 

r(a) dt 



ds(t — s) a x s 2 exp [ — — ) exp 



s 

T 2 



'dl£ 
dz 



A(s) (65) 



We denote 



B = C a 7 



<9iT 

<9z 



1 /&fT 



(66) 
(67) 



and make the Laplace transform of ( 1651) . Denoting the Laplace transform of 
the function A(t) as f3{p) 



A(t) <- /3(p) 



(68) 



we obtain after standard manipulations 



(p + c)/9(p + c) - A(0) = -5p 



!_ Q rf 2 /5(p + c) 



dp 2 



(69) 



Introducing a new complex variable 



z = p + c 



(70) 



we obtain an ordinary differential equation 



B(z - cf- a ^y- + = A(0) 



(71) 
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The problem we encouter with is to solve the homogeneous part of the equation 
( ffB . Introducing a new function u(z) with the help of the relationship 

d(3{z) 



dz 



(3(z)u(z) (72) 



we can cast the homogeneous part of the equation (ITTi) into an equivalent form 

du(z) 2 z , , 

— ; Vu{z) = -— - — (73) 

dz v ' B(z-cY- a v ' 



The latter is a Riccati equation. As is well known there is no general way 
to obtain analytical solution of the Riccati equation and the author of this 
chapter failed to obtain the solution for our particular case. In our opinion 
the approach encounters severe mathematical difficulties. We conclude that 
for the case of anomalous diffusion in complex media the NMR relaxation 
method has advantage over the PFG NMR because its mathemitical scheme is 
analytically tractable as the previous Sec. shows. The latter conclusion should 
be reconsidered because since this Chapter was submitted the problem of 
fractional Torrey equation has been solved in |31j . In the next Sec. we discuss 
some difficulties of physical character that still remain and prevent the NMR 
relaxation method to a become flawlessly accomplished one. 



6 Conclusion 

As is well known the main problem of the NMR relaxation method is that "the 
relaxation mechanism of the probe species needs to be known, and it is required 
that the intermolecular contributions to the relaxation can be separated from 
the intramolecular contributions" (see [27J and refs. therein). As for the first 
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part of the above argument the dipole-dipole relaxation mechanism dominates 
in very many cases and can be reliably supported upon. However the second 
part is really a problem. To put it in other words there are contributions to 
the spin-lattice relaxation rate constant from both translational and rotational 
diffusion of the tracer molecule 

(1/H) = + (l/Ti) rot (74) 



The input information for the approach developed in this chapter is (1/Ti) trans 
and the problem is how to separate it reliably from the rotationl contribution, 
i.e., to extract it from the experimentally observable value (1/Ti). For the case 
of ordinary diffusion (both rotational and translational) the problem is solved 
somehow because we know an explicit expression for (1/Ti) rat [17]. However 
in a complex media both translational and rotational contributions generally 
become some functions of the parameter a characterising the extent of inho- 
mogeneity. While the dependence (1/Ti) trans (a>) can be calculated explicitly 
(see [16J and Sec. 3 in this chapter) the dependence (l/Ti) rot (a) remains un- 
known. To develop mathematically the approach for calculation of the latter 
dependence along the line of [T7j basing upon the results of the paper [32] is 
quite feasible. However a physical problem still remains: whether the parame- 
ter a figuring in the fractional derivative of the equation for the translational 
diffusion should be the same as that for the rotational one? If they differ from 
each other what is the relationship between them? The author of this chapter 
does not know the answers to these questions at present. Untill this problem is 
resolved one should resort to some intuitive physical arguments. In our opin- 
ion it seems reasonable to assume that the rotational diffusion depends on the 
extent of inhomogeneity in a much narrow range than the translational one. 
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For instance in a porous media the translational diffusion is very sensitive to 
the radius of the pores in the whole range of this parameter while the rota- 
tional diffusion seems to become sensitive to it only when the radius of the 
pores becomes commensurable with the radius of the tracer molecule. If it is 
really so then in a wide range of the extent of inhomogeneity a c < a < 1 one 
can assume that (l/T 1 ) ro4 is independent of a and use the expression for it 
from [T7] obtained for ordinary rotational diffusion. This trick enanles one to 
extract the required value (1/T\) trans (a) from the experimentally observable 
one (1/31) (a) 

(l/r 1 ) traM (a) = (1/31) (a) - (l/31) ro t (75) 



The result of such procedure is the input information for the approach de- 
veloped in this chapter. If it is obtained then the below described strategy is 
feasible. 

The formula (161]) relates the FDC with the function H(x,a) fl39|) . The latter 
takes into account the contribution to the spin-lattice relaxation rate constant 
by anomalous translational diffusion and can in principle be extracted from the 
experiment as is described above. Besides a it depends on the dimentionless 
parameter x = lol,t\- In practice it is rather problematic to vary the Larmor 
frequency lu^. On the other hand the characteristic time (JB]) t\ = X2 ^^P de- 
pends on temperature and can be easily varied during the experiment. This 
fact enables one to retrieve the dependence of the contribution to the spin- 
lattice relaxation time due to translational diffusion (T\) trans [x, a) on x at 
different a from the experimental measurements. Approximation of the data 
obtained by a suitable analytical function of two variables yields H(x, a) and 
substitution of the latter into ( 16T1) yields the required dependence of the FDC 
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on the parameter a. 

We conclude that the formula (l6Tj) is an ingredient in solving the problem 
- how to retrieve the fractional the FDC from NMR relaxation data? The 
mathematical scheme for NMR relaxation method is analytically tractable in 
contrast to that of PFG NMR. In the latter case the attempts to develop a 
theory suitable for the description of the anomalous diffusion encouter severe 
mathematical difficulties. The result obtained in this chapter promotes the 
NMR relaxation method to become a powerful tool in solving the problem 
of experimental measuring the fractional diffusion coefficient. However much 
work still remains to be done for reliable separating the contribution of trans- 
lational diffusion into the spin-lattice relaxation rate constant from that of 
rotational diffusion in the experimentally measured values. This preliminary 
step is necessary for obtaining input information for the approach developed 
in this chapter. 
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